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Abstract 

The annihilation-creation operators a" are denned as the positive/negative fre- 
quency parts of the exact Heisenberg operator solution for the 'sinusoidal coordinate'. 
Thus are hermitian conjugate to each other and the relative weights of various 
terms in them are solely determined by the energy spectrum. This unified method ap- 
plies to most of the solvable quantum mechanics of single degree of freedom including 
those belonging to the 'discrete' quantum mechanics. 



1 Introduction 

The annihilation and creation operators are probably the most basic and important tools in 
quantum mechanics. Modern quantum physics is almost unthinkable without them. About 
eighty years after its birth, the list of exactly solvable systems in quantum mechanics |T] is 
quite long now, including those of the so-called 'discrete' quantum mechanics 0|. One 
natural question is that if these exactly solvable quantum mechanical systems also possess the 
algebraic solution method embodied in the annihilation-creation operators. We will answer 
the question in the affirmative and give a unified dynamical derivation of the annihilation- 
creation operators for most of the solvable quantum mechanics of single degree of freedom 
including those belonging to the 'discrete' quantum mechanics. 



The method is quite simple and elementary. One identifies a special function i](x) of the 
space-coordinate x, which undergoes 'sinusoidal motion' at the classical and quantum levels 
()2.10|) . (|2.14p . The latter is simply the exact Heisenberg operator solution for i](x). The 
function rj is the argument of the orthogonal polynomial ()2.3|) constituting the eigenfunc- 
tions of the system. The positive/negative frequency parts of the exact Heisenberg operator 
solution of the 'sinusoidal coordinate' rj(x) (j2.28J) . ()2.29|) give the annihilation-creation op- 
erators. This is essentially the same recipe as used by Heisenberg for solving the harmonic 
oscillator in matrix mechanics. To the best of our knowledge, the 'sinusoidal coordinate' was 
first introduced in a rather broad sense for general (not necessarily solvable) potentials as a 
useful means for coherent state research by Nieto and Simmons [3]. 

By similarity transformation in terms of the ground state wavefunction, the results of the 
present paper will be translated to those of the corresponding orthogonal polynomials. In 
particular, the exact Heisenberg operator solution of the 'sinusoidal coordinate' corresponds 
to the so-called structure relation [5J for the orthogonal polynomials. Our method provides 
the unified derivation of the structure relations for the Askey- Wilson, Wilson, continuous 
dual Hahn, continuous Hahn and Meixner-Pollaczek polynomials [S| based on the Hamilto- 
nian principle. These polynomials are the eigenf unctions of the various 'discrete' quantum 
mechanical systems [2J El and they are the deformations of the Jacobi, Laguerre and Hermite 
polynomials [7j. 

This paper is organised as follows. The general theory of the annihilation-creation op- 
erators is explained with typical examples; one from the ordinary quantum mechanics and 
two from the 'discrete' version in section two. Further explicit results are presented in sec- 
tion three. They include the (symmetric) Poshl- Teller, symmetric Rosen-Morse, Morse and 
x 2 + 1/x 2 potentials on top of the five examples belonging to the 'discrete' quantum mechan- 
ics mentioned above. Section four is for a summary and comments. In Appendix A, the 
necessary and sufficient condition for the existence of the 'sinusoidal coordinate' is analysed 
within the context of ordinary quantum mechanics. It turns out that those potentials having 
the 'sinusoidal coordinate' are all shape invariant [8J. In Appendix B, interpretation of the 
annihilation-creation operators within the framework of shape invariance is given. It is the 
mechanism underlying the solvability of all systems considered in this paper. Appendix C 
gives various definitions of the orthogonal polynomials, hypergeometric functions and their 
g-analogues E] . 
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2 General theory with typical examples 



The purpose of the present paper is to present a unified dynamical theory of annihilation- 
creation operators. It is applicable to most of the exactly solvable quantum mechanical 
systems of one degree of freedom, including the so-called 'discrete' quantum mechanics which 
are certain deformation of the solvable quantum mechanics j2j Oj. They satisfy certain 
difference equations instead of the second order differential equations. Generalisation to 
the systems of many degrees of freedom will be discussed elsewhere. The restriction to 
the solvable quantum systems is rather trivial and inevitable, since a system is obviously 
solvable if it possesses explicitly defined annihilation and creation operators and any one 
single eigenstate to work on. Then the entire set of exact eigenstates are easily and concretely 
generated. 

Except for the simple harmonic oscillator, which gives probably the only so far universally 
accepted example of the annihilation-creation operators, there are quite a wide variety of 
proposed annihilation and creation operators in the literature jHj • Historically most of these 
annihilation-creation operators are connected to the so-called algebraic theory of coherent 
states, which are usually defined as eigenstates of annihilation operators. Therefore, for a 
given potential or a quantum Hamiltonian, there could be as many coherent states as the 
definitions of the annihilation operators. 

Our new unified definition of the annihilation-creation operators is, on the contrary, 
based on the dynamical properties of a special coordinate, the 'sinusoidal coordinate' shared 
by a class of solvable dynamical systems discussed in this paper. A quantum mechanical 
system with a self-adjoint Hamiltonian TC is solvable (or solved) if the entire set of its energy 
eigenvalues {S n } and the corresponding eigenvectors {<fi n }, n = 0, 1, . . . are known: 

H(j) n = S n (f> n , n = 0, 1,...,. (2.1) 

As is well-known a quantum Hamiltonian (together with its 'discrete' analogue) has in general 
discrete as well as continuous spectrum. In this paper we will concentrate on the discrete 
energy levels only, either finite or infinite in number. Then, because of the one-dimensionality, 
the eigenvalues are not degenerate 

So < Si < ■ ■ ■ , (2.2) 
and the eigenvectors have finite norms ||0 n || 2 = N% < oo. (When normalised vectors are 
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needed we denote them by adding a hat, {0 n = f <p n /N n }, \\4> n \\ = 1.) This is the solution 
in the Schrodinger picture and the eigenvectors {<f) n } are usually expressed as functions 
{(p n (x)} of the space-coordinate x. For the majority of the solvable quantum systems, the 
n-th eigenfunction has the following general structure pQ: 

<f> n {x) = (f) (x)Pn(v(x)) (2.3) 

in which (po(x) is the ground state wavefunction. It has no nodes and we may choose it to 
be always real and positive. The second factor P n (rj(x)) is a polynomial of degree n in a 
real variable rj. We also take P n {i]) as rea l an d use a convention P_i(^) = 0. Reflecting 
the orthogonality theorem of the eigenvectors of a self-adjoint Hamiltonian, {P n (r])} form 
orthogonal polynomials with respect to a weight function (measure) 

(f>o(x) 2 dx oc w(r))dr). (2.4) 

Throughout this paper we follow the definition and notation of Szego's book 7\ for the 
classical orthogonal polynomials and the review by Koekoek and Swarttouw [0] for the Askey- 
scheme of hypergeometric orthogonal polynomials and its g-analogue, which are deformations 
of the classical orthogonal polynomials. 

There are certain exactly solvable quantum systems, for example, the one- dimensional 
Kepler problem, etc, for which the general form of eigenfunction (|2.3j) does not hold. For 
them, the present unified theory does not apply. See Appendix A. 

Our main claim is that this r](x) undergoes a 'sinusoidal motion 1 under the given Hamil- 
tonian 7i, at the classical as well as quantum level, by mimicking the simple harmonic oscilla- 
tor. This fact is the basis of our dynamical and unified definition of the annihilation-creation 
operators. To be more specific, at the classical level we have 

{H, {H,n}} = -vMK) - R-i{H) (2.5) 

in which the canonical Poisson bracket relations are defined for the canonical coordinate x, 
its conjugate momentum p and for any functions A(x,p) and B(x,p) as: 

f) A f)R f) A f)R 

{«,,} = !, {*,*} = {P , P } = 0, {A,B} = (2.6) 

The two coefficients Rq and R-\ are, in general, polynomials in the Hamiltonian 7i. The 
effect of R-i is to shift the origin of rj(x) by a quantity (possibly) depending on Tt. It is 
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convenient to introduce a shifted sinusoidal coordinate rj(x) 

fj( x ) d ^ f v (x) + R^(H)/R (H) (=> {ft, {ft, 77}} = -fj R (H)) . (2.7) 
The relation (|2.5jl would allow to evaluate the multiple Poisson brackets of 77 with ft easily: 
ad ft 77 = {ft, 77}, (adft) 2 77 = {ft, {ft, 77}}, (adft) n 7/ = {ft, (adft)™^}, (2.8) 
which leads to a simple sinusoidal time-evolution: 

z — ' n! 

n=0 

= -{ft, r/} ^^p] + 77(^)0 cos [t ] • (2-9) 
V tto\"o) 

In the original variable it reads 

sin \tJR (H )~ 

V ( X ; t) = -{ft, V } l )-^- n ~ R-lW)/RoW) 

+ (n(x) + R^(n )/Ro(^o)) cos [V^o(fto) ] , (2.10) 

in which 77(2)0 ( ^(^)o) an d {ft, t/}o are the initial values (at t — 0) of these variables and 
fto denotes the value of the Hamiltonian (the energy) for these initial data. In general, 
the frequency of the simple oscillation a/ -Ro(fto) can depend on the initial data. This is 
the reason why we call 77(0;) the 'sinusoidal coordinate' avoiding the more appealing but 
misleading "harmonic coordinate". 

At the quantum level with the canonical commutation relations 

[x,p]=i, [x,x] = \p,p] = 0, H = l, (2.11) 

the formula corresponding to (|2.5j) reads 

[ft, [ft, 77] ] =r ] Ro(H) + [ft, 77] i?!(ft) + itU(ft). (2.12) 

In other words, the multiple commutators of ft with 77 form a closed algebra at level two. 
Here, the quantum coefficients Rq and R^i could differ from the classical ones by quantum 
corrections. But we use the same symbols since there is no risk of confusion. Obviously 
i?i(ft) is the quantum effect. As in the classical case, the multiple commutators of 77 with 
ft 

adft77= [ft,77], (adft) 2 77 = [ft, [ft,7/]], (adft) n 77 = [ft, (adft)"" 1 ^], (2.13) 
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can be easily evaluated from (j2.12|) . This leads to the exact operator solution in the Heisen- 
berg picture: 

e un V (x)e- un = jr (adH) n V , 

n=0 U ' 

= |H -' wl MJFMH)^- l(H)WI 

in which the two "frequencies" a±(TL) are 



ol±{H) = {Ri{H) ± ^Ri(H) 2 + 4R (H))/2, (2.15) 
a+(H) + or_(W) = Ri(H), a + (H)a^(H) = -R (H). (2.16) 

If the quantum effects are neglected, i.e. R\ = and H — > 7i (in the r.h.s.), we have 
a + = — ct_ = a/ Ro(Hq), the above Heisenberg operator solution reduces to the classical one 
(J2.1(J|) in terms of the quantum-classical correspondence: 

[A, B]/ih — > {A, B}, (/wO). (2.17) 

The above exact operator solution looks slightly simpler if the shifted sinusoidal coordinate 
is used 

itn~, \ -an ^ ~, ^ M+(W)t - e^ (w)i . -a_(ft)e ia +W* + a+(H)e iQ -( w )' , oio , 
e-,(x)e a+(?0 + « + (70-«.(7fl " ^ 

Like the exact classical solution ()2.1()jl . the exact quantum solution (|2.14jl contains all the 

dynamical information of the quantum system. One can, for example, determine the entire 

discrete spectrum {S n } by following Heisenberg and Pauli's arguments for the harmonic 

oscillator and the Hydrogen atom. Let us first note that the ground state energy £ = 

is known explicitly, because of our choice of the factorised form of the exactly solvable 

Hamiltonian (see examples below): 

H = A*A/2, ^0 O = O^H0 O = O, £ = 0. (2.19) 

Let us apply (|2.14|) to the n-th eigenvector (p n : 

e it(n - £n) ri(x) <f> n 

'[HMz)]<t>n + {-V(x)a-(£ n ) + J R_ 1 (^,)/a + (^))0„)e ia +^)7(a + (^,) - e*_(£ n )) 
+ (-[H,v(?)]<i>n + (v(x)a + (S n ) - J R„ 1 (^,)/«-(^))0n)e ia - (£ " ) 7(a + (^) - a_(£ n )) 
- (R^{E n )/R Q {E n ))^ n . (2.20) 
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Since the r.h.s. has only two different time dependence except for the constant term, the 
l.h.s. can only have two non-vanishing matrix elements when sandwiched by <f) m , except for 
the obvious <p n corresponding to the constant term. In accordance with the general structure 
of the eigenfunctions (|2.3jl . they are 4> n ±i- 

(<Pm\v(x)\<Pn) = 0, for m^n±l,n. (2.21) 

This imposes the following conditions on the energy eigenvalues 

£n+l - £n = a +(£n), £n-l ~ £ n = Ot-(£ n ). (2.22) 

Likewise we obtain the 'hermitian conjugate' conditions 

£n — £n~i — oi+(£ n -i), £ n — £ n +i = ot-(£ n+ i) (2.23) 

relating these three neighbouring eigenvalues. These overdetermined conditions ()2.22)1 . (J2.23)) 
and £q = determine the entire energy spectrum {£ n } completely for each Hamiltonian. The 
consistency of the procedure requires that the second term on r.h.s. of ()2.2()j) should vanish 
when applied to the ground state </> : 

-[H,v(x)]<l>o + (v(x)a + (0) - itU(O)/a_(O))0 o = 0, (2.24) 

which could be interpreted as the equation determining the ground state eigenvector O m 
the Heisenberg picture. If the number of the discrete levels is finite (M + 1) a corresponding 
condition must be met that the first term on r.h.s. of ()2.2()j) should not belong to the Hilbert 
space of normalisable vectors when applied to the highest discrete level eigenvector <pM'- 

\\[H,r](x)}(f)M + (-T)(x)a-(£ M ) + R-i(£m)/ol+{£m))4>m\ \ = oo. (2.25) 

As is clear by now, ()2.14j) and 1)2.21)1 are the physical embodiment of the three term 
recursion relations satisfied by any orthogonal polynomial of single variable: 

r)P n {ri) = A n P n+1 ( V ) + B n P n (rj) + C n P n ^(rj). (2.26) 

The coefficients A n , B n and C n are also real for a real polynomial P n {v)- F° r t ne systems 
treated in this paper, that is those having the general structure of the eigenvectors ()2.3)) . 
this implies the three term recurrence relations of the eigenfunctions 

r)(x)4> n (x) = A n (f) n+1 (x) + B n (j) n (x) + C n 0„_i(x). (2.27) 
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At the same time the above arguments and the treatment of the harmonic oscillator by 
Heisenberg clearly show that the operator coefficient of e Jto -( w ) on the r.h.s. of the Heisenberg 
operator solution (J2.14J1 is the annihilation operator, that is, acting on n it produces a 
state (j) n -i- Likewise, the operator coefficient of e 4to +^^ of the Heisenberg operator solution 
()2.14|) is the creation operator. Thus we arrive at a dynamical and unified definition of the 
annihilation-creation operators: 

e m r](x)e- ltn = a M (H, r])e ia+{H)t + a { -\n,v)e ia - {n)t - R-i(H) / R (H) , (2.28) 

= (±[H,v(x)] =f (tj(x) + R^(H)/Ro(H))a T (H)} / {a + (H) - ou.{H)). (2.29) 
When acting on the eigenvector (f> n , they read 

« (± Vn(x) = . ±1 p ([H, V (x)} + (£ n - £ nTl ) V (x) + *- l{£n l )<p n {x). (2.30) 

Before going to the detailed discussion of the annihilation-creation operators for various 
Hamiltonians in section |2~T1 and sectional let us analyse annihilation-creation operators in 
a more general context. A minimal requirement for annihilation-creation operators is the 
following: 

(0) : annihilation-creation operators map (fi n to 4> n -i and 4> n+ i (up to an overall 
constant), respectively. 

It should be stressed that there is no a priori principle for fixing the normalisation of the 
operators. Sometimes it is convenient to introduce the annihilation-creation operators with 
a different normalisation 

a /(±) drf a (±)( a+ ( W ) - a _(ft)) = ±[H,r](x)} =F (r]{x) + R^{n)/R (n))a T (H), (2.31) 

which gives 

a^VnO*) = ±([H,v(x)} + (£ n - £ nTl ) V (x) + ^ £n ] )<j> n {x). (2.32) 

In the r.h.s. the coefficients of the operator r](x) and the identity operator depend on n in 
general. 
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The annihilation-creation operators of the harmonic oscillator have several remarkable 
properties: 

(i) : Annihilation/creation operator is the positive /negative frequency part of 
the Heisenberg operator of the 'sinusoidal coordinate'. 

(ii) : (annihilation operator)^=(creation operator). 

(iii) : H = const. x (creation operator) (annihilation operator). 

The first property (i) is the principle leading to our unified definition of the annihilation- 
creation operators as in ()2.28|) - (J2.29|) . Next we show that they are hermitian conjugate to 
each other. That is, they satisfy the property (ii), too. By using the three-term recursion 
relation, we obtain 

e l 'V~ im 0n = e^+^A^+x + B n cj) n + e^-^OA-i. (2.33) 
Comparing this with (|2.28|) and (|2.29|) . we arrive at 

a<+Vn = A n( p n+1 , a^<j) n = C n n _ x , R„ 1 (S n )/R (£ n ) = -B n . (2.34) 
Here use is made of the facts 

ot+(S n ) = £ n +i — £n, ot-(£ n ) = £ n -\ — £ n , 

R\{£n) = £n+l + £n-l — 2£n, Ro(£n) = ~ (£n+l ~ £n)(£n-l ~ £n) , (2.35) 

and that ot±(H) and Ri(H) (i = 1,0, —1) are hermitian. Hermitian conjugate of is 

a Ht = (a + (H) - a^{H)Y l ([H, rj(x)] + a + (H) (rj(x) + R^(H) / Ro(H))\ (2.36) 

and its action on <p n is 

><->U n = (a + (H) - u_{H)Y l (Ur](f> n - r]£ n <f> n + a + (W)(# n + R^(£ n ) / Ro(£ n )(j> n ) 
.(H) -a_{H)Y l ({£ n+l -£ n )A n <f) n+1 



(a+{H) - a -(H)) (£ n +2 ~ £ n )A n <j) n+ i 

A n n+1 = a^4> n . (2.37) 
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Therefore are hermitian conjugate to each other, aHt = a (+). 

This also means that 

e - ito _(H) a (+) = a (+) e ito + (W) j e -ita + (W) a (-) = Q (-) e to_(W) ) (2>38) 

reflecting the obvious hermiticity of the l.h.s. of (|2.28j) . Note that a ,<y± ^ are not hermitian 
conjugate to each other, a'(~" 7^ a ,<y+ \ in general. 

In the special case of the equi-spaced spectrum £ n = an (a : constant), to which many 
interesting examples belong including the harmonic oscillator and its deformation, we have 
a±{H) = ±a, Ri(H) = 0, Rq(H) = a 2 and 

2aa (±) = a /(±) = ±[H, r}{x)\ + ar}{x) + R^(n)/a. (2.39) 

For the simplest harmonic oscillator 7i = (p + ix)(p — ix)/2, we have n(x) = x, Rq = 1, 
Ri = = and [H, x] = —ip and 2a^ + - ) = x — ip, 2a^ = x + ip, which differ from the 
conventional ones by a factor y/2. 

In contrast to the above two properties, the third the property (iii) of the annihilation- 
creation operators is achieved by a very specific modification of the definition as follows: 

a "(-) d£f a (_) a „(+) dcf a , /( _ )t = /(W)a (+) > (2 40) 

where /(W) = j^H)^ is an as yet unspecified function of TC. Then we have 

a"(+V'(-Vn = f(£ n ) 2 A n ^C n 4> n . (2.41) 

If /(W) is chosen to satisfy 

f{£ n f = £ n /(A n ^C n ) ( = g{n)), (2.42) 

then we obtain 

n = a "i+) a "(-), (2.43) 



Such an operator f(7i) function can be constructed as f(7i) = ygiN), where the number 
(or level) operator M is defined by 

N<t> n = n<j) n . (2.44) 
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This operator Af can be expressed in terms of the Hamiltonian, for example, 

£ n = an J\f=H/a, (2.45) 

£ n = an 2 + bn A/" = (y/AaH + b 2 - 6) /(2a), (2.46) 

£ n = a(q' n -\){l-bq n ) q M = (H/a + b + 1 - y/(H/a + b + l) 2 - 46)/(26), (2.47) 

where a,b,q are constant (6 > in the second equation, 6 < 1 in the third equation). This 
a"^ satisfies the property (ii) by definition, but the property (i) becomes ugly or unnatural 
in general. 

Let us close the general theory by a brief discussion of the coherent states. One definition 
of the coherent state if) is the eigenvector of the annihilation operator (AOCS, Annihilation 
Operator Coherent State): 

a (~V = \if), AG C. (2.48) 

In terms of the simple parametrisation if) = Y^=o c n<t>n{.x) (with c = 1 as normalization) 
and the formula (|2.34j) . we arrive at 

if) = -0(A, X) = o (z) rrn r ■ Pn(jl{x))- (2-49) 
n=0 llfc=l°n 

For the equi-spaced spectrum £ n = an (a : constant), the coherent state has the property of 
temporal stability 

e itn if){X,x) = if){e iat X,x). (2.50) 

It should be remarked that the concrete form of the AOCS depends on the specific normal- 
isation of the annihilation operator. For the annihilation operators a'(~\ a"^~' and others, 
we denote the corresponding coherent states as ip', if>", etc. Which coherent state is useful 
depends on the physics of the system. 

2.1 Some typical examples 

Now let us look at typical examples |H El to show the actual content of our new unified 
theory of annihilation-creation operators. In their pioneering work, Nieto and Simmons [3] 
treated four solvable cases, those discussed in sections 12.1.11 13.1.11 13.1.31 and 13. 1.41 Some of 
our results were reported in jl]. Here and throughout this paper we put the dimensionfull 
quantities as unity, including the Planck's constant. 
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2.1.1 l/sin 2 x potential, or symmetric Poschl-Teller potential 

The first example has the 1 / sin 2 x potential, which is the one-body case of the well-known 
Sutherland model [TO"! lllj. This provides the simplest example of the annihilation-creation 
operators depending on n. The corresponding coherent state (|2.72j) or (|2.75jl had not yet 
been known, to the best of our knowledge. The system is confined in a finite interval, 
say (0, 7r) and it has an infinite number of discrete eigenstates. Although this potential is 
a special (g = h) case of the Poschl-Teller potential discussed in section 13.1.21 it merits 
separate analysis. The Hamiltonian, the eigenvalues and the eigenfunctions are as follows: 

H = (p - ig cot x)(p + ig cot x)/2, ( ^ 2H + g 2 = p 2 + g{g - 1)/ sin 2 x) , (2.51) 
£ n = n{n/2 + g), n — 0,1,2,..., g > 0, < x < n, rj{x)= cosx, (2.52) 
(j) n (x) = {sin x) 9 Pjf'® {cosx), (3 = g- 1/2, (2.53) 

in which P„ ^ (rj) is the Jacobi polynomial (jC8|) and P^ 3 '^ {rj) is proportional to the Gegen- 
bauer polynomial d' 8+1//2 ' ) {rj) (jC.9|) 

(P + 1)„ (20 + 1)„ ' [ ' 1 

Hereafter we often use the Pochhammer symbol {a) n , see (jC.lJl . 
It is straightforward to evaluate the Poisson brackets 

{H,cosx} =psmx, {H,{H, cosx}} = -cosx 2H', H' = H + g 2 /2, (2.55) 

leading to the solution of the initial value problem: 

r / 1 sinrtv / 2HT 

cosx(t) = cosx 

(0)cos[tV 2 ^oJ -p(0)sinz(0) — 1 f UJ . (2.56) 

/2H.Q 



It is straightforward to verify | cosx{t)\ < 1. The corresponding quantum expressions are 

[H, cosx] = i sin xp + cos x/2, (2-57) 
[H,[H, cos ar] ] = cos x{2H' - 1/4) + [H, cos x] , (2.58) 



a ± {H) = 1/2±\W. (2.59) 
The exact operator solution reads 

e im cosxe~ %m = (i sin x p + cos x/2) ; 

2 v / 27f 

-aJH\e^W + amy™-™ , N 

+ cos x — ; + v ' . 2.60 

2V2H' 
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The annihilation and creation operators are 

a'(±) =a (±) 2v / 27f = ±isinxp + cosxv / 27f = ±sinx-^- + cosxv / 27^. (2.61) 

dx 

It is now obvious that they (a'^) are not hermitian conjugate to each other. The square 
root sign is neatly removed when applied to the eigenvector <p n as 2S n + g 2 = (n + g) 2 : 

a '^Vn = _ sin x— ^ + (n + g) cos x (j) n = {n + P)(p n _i, (2.62) 
dx 

a' i+) 4> n = sinx^ 1 + (n + g) cosx0 w = ^ " +^ ^2g + \ 9 ^ n+l ' ( 2-63 ) 

This is a rule rather than exception as expected from the relations between the neighbouring 
energy levels, ()2.22|) and (j2.23|) . The right hand sides are the results of the application. 
In particular, when acting on the ground state 4>o, the annihilation (a'^) and creation 
(a'( + )) operators are proportional to the factorisation operators A and At of the Hamiltonian 
7i = A^A/2, respectively: 

d x , / d 



a 



'(-), 



>o = (— sinx- — h g cosx)0 o = — sinx(- g cot x)0 o = V \x)A4>o = 0, 



dx dx 
a'( + Vo = (sins- — hgcosx)0 o = — sinx(— gcotx)</> — r}'(x)A^(j)o- (2.64) 

(X Jij KAjJLj 

In a rough sense, the factor n'(x) = — sinx, in the creation operator, compensates the 
downward shift of the parameter (g) caused by A^ . Similar situations are encountered in 
all the other quantum systems. In particular, for the systems with equi-spaced spectrum 
£ n = an (a : constant), the factorisation of a'* - -* and a'^ into A and A^ is n-independent, 
(|2.94|) . (|3.14j) . (|3.15j) . (|3.92j) . Their significance will be discussed in some detail for the 
'discrete' quantum mechanics cases in Appendix B. 

The following interesting commutation relations ensue from (|2.62j) and (J2.63|) : 



[H,a l{±) ] = ±(Vm>a'^ +a /(±) v / 27f)/2, (2.65) 
[a'(-),a' {+) ] =2v / 27f, (2.66) 
a 'H a '(+) + a'^a'^ =AH + 2g. (2.67) 

The relation ()2.67|) could be accepted as a substitute of the property (iii) of the annihilation- 
creation operators discussed in page M 

By the similarity transformation in terms of the ground state wavefunction <po{x) = 
(sinx) 3 , we obtain the so-called shift down and up operators for the Jacobi (Gegenbauer) 
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polynomial (/3 = g — 1/2): 

down : (l-r ] 2 )^-P^ l3 \r ] )+nr ] P^ l3 \r ] ) = (n + P)Pjf^\r ] ), (2.68) 
ar] 

up : -(1 - V 2 )^ m (v) + (n + 2g) V Pr\ V ) = P&fiv)- (2-69) 

As expected they are the same Jacobi polynomials of degree n — 1 and n + 1. It should 
be stressed that these shift down-up operators are naturally derived from our annihilation- 
creation operators without assuming the explicit form of the three term recursion relation. 

The coherent state ip (I2.48J) is 

i/,(x) = o (x) f;(2A) n ( ^ + 3 ^ )w P B ^(cosx) = o (x) f;(2A)»^±^CW(oosx), (2.70) 

n=0 ^ n=0 ^ ®> n 

where we have used a^~'(j) n = a'(~> (j> n / (2(n + g)), ()2.62|) and (|2.54|) . A generating function of 
the Gegenbauer polynomials (7: arbitrary) 



n=0 

gives a concise expression of the coherent state i/j 



-2A sin 2 x 



^) = (l-2Acosx)-- 1 2 F 1 f ( ^ + 1)/2 ^ /2 + 1 - — - ). (2.72) 
^ K J K J : V £ + 1/2 (l-2Acosx) 2 / V ; 

Here %Fi is the hypergeometric function ()C.3|) . For the annihilation operator a'(~\ the 
corresponding coherent state ip' is 



f(x) =M*)J2nrhv P n m ( C0SX ) = ^o(x)^-Vc^(cosx). (2.73) 

n=0 V 3 + ^ n=0 ^ 9)n 

A generating function of the Gegenbauer polynomials [6 , 

|^.(^^( 9+ - /2 |^), (2,4) 

gives a concise expression of the coherent state ip': 

i>'(x) =T(g + 1/2) e Xcosx (\/2) 1 / 2 - 9 V^J g - 1/2 (\smx), (2.75) 
in which J a (x) is the Bessel function (IC.5I) . 
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2.1.2 Deformed harmonic oscillator =>- Meixner-Pollaczek polynomial 

The deformed harmonic oscillator is a simplest example of shape invariant 'discrete' quantum 
mechanics. The Hamiltonian of 'discrete' quantum mechanics studied in this paper has the 
following form [2] (with some modification for the Askey- Wilson case in section l2.1.3j) : 

U = (^V{x)e p ^V(x)* + ^/V(x)*e- p ^V(x) - V(x) - V(x)*^j / 2. (2.76) 

The eigenvalue problem for TC, TC<fi = £<fi is a difference equation, instead of a second order 
differential equation. Let us define S±, T± and A by 

S + = e p / 2 v^F, SL = e~*/ 2 y/V(x), SX = ^V{x)e p l\ Si = ^/V{xfe^ 2 , (2.77) 
T + = SlS + = ^V{x)e p ^/V{xy-, T_ = SlS- = ^/V{xY e~ p ^/V{x) , (2.78) 
A = t{S + - S-), A ] = -i(Sl - Si). (2.79) 

Then the Hamiltonian is factorized 

H = (T+ + T_ - V(x) - V{x)*)/2 = (Si - Sl)(S+ - 5_)/2 = A* A/2. (2.80) 

The potential function V(x) of the deformed harmonic oscillator is 

V(x) = a + ix, — oo < x < oo, a > 0. (2-81) 

As shown in some detail in our previous paper [2], it has an equi-spaced spectrum and 
the corresponding eigenfunctions are a special case of the Meixner-Pollaczek polynomial 

p^(x^) CUB, 

S n = n, n = 0,1,2,..., (2.82) 

4>o(x) = \/r(a + ix)T(a — ix), rj(x) = x, (2.83) 

<P n (x) = <h(x)P n (x), P n (x) = Pi a \x; f ), (2.84) 

which could be considered as a deformation of the Hermite polynomial. 
The Poisson bracket relations are 

{H,x} = -Va 2 + x 2 sinhp, {H, {H, x}} = -x, (2.85) 

leading to the harmonic oscillation, 

x(t) = x(0) cost + ^a 2 + x 2 (0) sinhp(O) sint, (2.86) 
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which endorses the naming of the deformed harmonic oscillator. The corresponding quantum 
expressions are also simple: 

[H, x] = -i(T + - T_)/2, [H, [H, x]]=x, (2.87) 
e m xe~ im = x cost + i[H,x] sint = x cost + (T + - T_)/2 sin t. (2.88) 

The annihilation and creation operators are 

a'W = 2a {±) =x±{H,x}=xt i{T+ - T_)/2, (2.89) 

which are hermitian conjugate to each other. These operators were also introduced by De- 
gasperis and Ruijsenaars ^2] by a different reasoning from ours. By similarity transformation 



in terms of the ground state wavefunction 4>q(x) = yT(a + ix)T(a — ix), we obtain 

<p Q (x)- 1 a' (±) o (x) = x =f i(V{x)e p - V{x)* e - p )/2 . (2.90) 

The action of the annihilation creation operators on the eigenvectors 

a '(-)0 n = (n + 2a - l)0 n _ 1; a' (+) 0™ = (n + l)0„+i (2.91) 

is consistent with the three term recurrence relation of the Meixner-Pollaczek polynomial: 

(n + ; f ) - 2xP(«) ( x; + + 2a- l)P^(x ; f ) = 0. (2.92) 

From these it is easy to verify the su(l, 1) commutation relations including the Hamiltonian 
H: 

[H, a' (±) ] = ± a ,(±) , [a ,( - } , a /(+) ] = 2(W + a). (2.93) 
It is interesting to note that are factorised by the factors of the Hamiltonian Ti = A^A/2, 

4a ( ~) = XU, 4a (+) = A*X, (2.94) 
X = S + + S-, X^ = Sl + SL (2.95) 

These X and compensate the shift of the parameter a caused by A^ and A. See Appendix 
B for more details. 

The coherent state (|2.48|) . (|2.49|) . is simply obtained from the formula (|2.91|) and a'^~> = 
2aH : 

^(x)=Mx)f^^Pi a) (x;l). (2.96) 

n=0 
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A generating function of the Meixner-Pollaczek polynomial 



(2a 

n=0 v 

gives a concise expression of the coherent state if>: 



${x) = M*) <? iX i F i C o %X ) > ( 2 - 98 ) 



2a 

in which 1F1 is the hypergeometric function ()C.3j) . 
2.1.3 Askey- Wilson polynomial 

The Askey- Wilson polynomial belongs to the so-called g-scheme of hypergeometric polyno- 
mials 6j. It has four parameters 01,02,03,04 on top of q (0 < q < 1), and is considered as 
a three-parameter deformation of the Jacobi polynomial. As shown in our previous papers 
[312], it also describes the equilibrium positions of the trigonometric Ruijsenaars-Schneider 
systems based on the BC root system 13J. Thus, as a dynamical system, it could be called 
a deformed Poschl- Teller potential or one body case of the trigonometric BC Ruijsenaars- 
Schneider systems. The quantum-classical correspondence has some more subtlety than 
the other 'discrete' quantum mechanical systems treated in section l3~2*l because of another 
'classical' limit q — > 1. 

The factorised Hamiltonian of the Askey- Wilson polynomial has a bit different form from 
that of the Meixner-Pollaczek polynomial (j2.76l) : 

H = (^V{z)q D ^V{zY + VvWq- D ^/V(zj-V(z) - V{zf) I 2, (2.99) 
with a potential function Viz): 

t r / x rL=l(l ~ a j Z ) ix n r> def d .d , nlnn > 

v(t > = (i-^d -,»') ■ z = e • 0<x<7r ' D = z Tz = -'T x =p - (2 ' 100) 

We assume —1 < 01,02,03,04 < 1 and 01020304 < q. This Hamiltonian is also factorised 
7i = A^A/2, where A and A^ are given in ()2.77j) - (j2.79|) with the replacement V(x) =>- V(z), 
e ±p/ 2 q ±D / 2 ^ e tc. The eigenvalues and eigenfunctions are [210]: 

£ n = (g- n -l)(l-oio 2 o 3 o 4 g n - 1 )/2, n = 0,1,2,..., (2.101) 



",.(•'•>- J n4 (Z / ?)M , r^T^K ' V{x) = t±f^ = ™.r. (2.102) 



4> n (x) = 4> (x)P n (cosx), P n (rj) = f Pn(^;ai,a 2 ,o 3 ,a 4 |g), (2.103) 
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in which p n (r) ; a±, 0,2, 03, a^\q) is the Askey- Wilson polynomial ()C.14|) . 

The presence of the g-factor has only superficial effects at the classical level with the 
Hamiltonian (7 = logg): 

H c = y/V c (z)V c (z)* coshyp - (V c {z) + V c (z)*)/2, V c (z) = ^ J ^~ 2 y Z \ (2-104) 



{7i c ,cosx} = 7y rij = i(l _ a 3 z ) llj=i( 1 _ a j/ z ) I (4sinx) sinli7p, (2.105) 

{H c , {H c , cos x}} = - cos x R (H C ) - R^iHc), (2.106) 

Ro(H c ) = 7 2 (K 2 + CiH c + c 2 ), R-i{H c ) = -i 2 (c 3 H c + c 4 ), (2.107) 

with coefficients Ci,..,c 4 : 

ci = l + 6 4 , c 2 = (1 - 6 4 ) 2 /4, c 3 = (61 + 6 3 )/4, C4 = (l-6 4 )(6i-63)/8. (2.108) 
Here we use the abbreviation 



4 

b t = ^ dj, 63 = ^ a j a kO>i, h = (2.109) 

l<j<4 l<j<fc</<4 j=l 

The corresponding quantum expressions are 

[H,cosx] = (q- 1 - l)(z -1 (l - qz 2 )T + + z(l - gz~ 2 )T_)/4, (2.110) 

[W, [W.cosrr]] = cosxR {H) + [W, cosar]fli(W) + R-i(H), (2-111) 

i?o(W) = ffGf 1 - 1) 2 (W 2 - (1 + g" 1 ) V*) , (2.112) 

R X (H) = qiq- 1 - l) 2 W, W = H + {1 + q-%)/2, (2.113) 

R^(H) = -q{q- 1 - l) 2 ((6i + q-%)H/A + (1 - q~%)(h - 6 3 )/s). (2.114) 

The two frequencies are: 

<*±(H) = (q- 1 - 1) ((1 - q)W ± (1 + q)V(H>) 2 -q- 1 h) / 2, (2.115) 

in which 

H'(j) n = {q~ n + hq n - l )/2 n , ((H') 2 - q-%)<P n = {q~ n - hq n ~ l ) 2 /A n . (2.116) 
The annihilation-creation operators are: 

a (±) = (±(g- x - l)(z _1 (l - qz 2 )T + + z(l - g^ 2 )T_)/4 

=F cosxa T (ft) ± i^iC^aiCW) -1 ) / (a+(W) - • (2.117) 
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Their effects on the eigenvectors are: 



° K = 2(1 - b^Wl - 6 4g ^-i) ^ (2 ' 118) 



a " = 2(l-6 4 ^-i)(l-6 4 ^) (2 ' 119) 
which are consistent with the three term recurrence relation of the Askey- Wilson polynomial. 
The 'annihilation-creation' operators on the polynomial P n (cosx) read 

0o(^) _1 4>o(x) ■ P n (cOSx) 

= ~ p (±(q- 1 ~ 1) (z-\l - qz 2 )V{z) q D + z(l - qz~ 2 )V \z)* q~ D ) / 4 

± (£, - S nT i) cosx ± ) P n (cos x) . (2.120) 

The coherent state is 

1/ \ At \P (2A) n (aia 2 a 3 a 4 ; g) 2n , , 
ip(x) = fo{x)}- r = p —P n {cosx). 2.121 

We are not aware if a concise summation formula exists or not. 

3 Various results 

In this section we will briefly present many interesting results on the annihilation-creation 
operators, their algebraic properties and coherent states, etc, for various exactly solvable 
quantum mechanical systems, including the 'discrete' quantum mechanical systems. All 
these solvable systems share shape invariance [El E 121 which is a purely quantum mechanical 
notion that guarantees quantum solvability. But that property plays no active role in the 
present theory. On the other hand, the 'sinusoidal motion' exists at the classical and quantum 
levels. In Appendix A we will show that a system realising the exact sinusoidal motion is 
quite limited and that all of them belong to the known shape invariant systems. In other 
words there are shape invariant systems that do not have sinusoidal motion. 

3.1 Ordinary quantum mechanical systems 
3.1.1 x 2 + 1/x 2 potential 

When a centrifugal barrier (a 1/x 2 potential) is added, the harmonic oscillator keeps its exact 
solvability, but the particle is restricted to a half line, either x > or x < 0. This is the one- 
body case of the well-known Calogero model [TUl HTJ . The eigenfunctions are described by 
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the Laguerre polynomial and the annihilation-creation operators within the su(l, 1) scheme 
are well-known |3J • Our unified theory predicts these operators naturally. The coherent 
and squeezed states in the su(l, 1) were already known [Tl]. Its Hamiltonian, the eigenvalues 
and the eigenfunctions are: 

TC = (p + ix — ig/x)(p — ix + ig/x)/2, < x < oo, g > 0, (3.1) 
g n = 2n, n = 0,1,2,..., 77(2) = x 2 , (3.2) 
0„(x) = n (x^) = e-^VLf (x 2 ), /5 d = f 0-1/2, (3.3) 

in which Ln(j]) is the Laguerre polynomial ()C.7|) . 
The Poisson bracket relations are simple 

{H,x 2 } = -2px, {H,{H,x 2 }} = -A(x 2 -H-g), (3.4) 

leading to the simple sinusoidal motion 

x 2 (t) =x 2 (0) cos2t + (l-cos2t){H + g) +p{0)x(0)sm2t. (3.5) 

It is straightforward to verify x 2 (t) > 0. The quantum theory is almost the same as the 
classical one: 

[U,x 2 ] = -i(xp + px), [H, [H, x 2 } } = A(x 2 - H'), W = H + g + 1/2, (3.6) 
e m x 2 e- itn = x 2 cos 2t + (1 - cos 2t)H' + (xp + px)/2 sin 2t, (3.7) 

which leads to the following annihilation and creation operators: 

a (±) = (x 2 -H')/2Tt{xp + px)/4= f(^ T x) 2 - 9 ^ 9 ~^ ) / A. (3.8) 

V ax x 2 / 

The action of these operators are (/3 = g — 1/2) 

a { - ] <pn = -(n + P)(f>n-i, a {+) <p n = -(n + l)<j> n+1 , (3.9) 
which are consistent with the three term recurrence relation of the Laguerre polynomial 

(n + l)L^ + M + (n-2n-fi- l)Lf{r,) + (n + flL&fo) = 0. (3.10) 
From these follow the su(l, 1) relations 

[^,a (±) ] =±2a {±) , [a(-\a i+) }=H' = H + g + 1/2. (3.11) 



20 



The coherent state (AOCS) ()2.49|) is obtained simply as 



^ = J2(JTWJ^ = ^ x ) e ^2 X ym U (W=A), (3.12) 



n=0 



in which a generating function of the Laguerre polynomial [H] 

oo n 

^(a + l)„ Va + 1 / 

n=0 v ' 

and ()C.5|) are used. 

The annihilation and creation operators (|3.8jl are factorised by the factors of the Hamil- 
tonian H = AU/2: 

4a ( ~> = + x + - V A, A = i(p-ix + ig/x), (3.14) 

4a (+) = A* ■ (—J- + x + -), = -z(p + ix - z#/x). (3.15) 

This is a degenerate case of (|3.92|) and the action of the other factors is to compensate the 
parameter shifts caused by the operators A and A: 

(~~(L + x + ~) < t ) ^ x '^9) = 2 <fin(x;g + 1), (3.16) 

(^-+x + ^)<j) n (x;g + l) = {2n + 2g + l)(t )n {x-g). (3.17) 
\dx xJ 

3.1.2 Poschl- Teller potential 

The Poschl- Teller potential has two parameters g and h and its eigenfunctions are related to 
the Jacobi polynomial. It is the one body case of the BC type Sutherland systems [TP] ITT]. 
Its Hamiltonian, the eigenvalues and the eigenfunctions are: 

7i == (p — igcotx + ih tan x)(p + igcotx — ihtanx)/2, 0<x<n/2, (3.18) 
£ n = 2n(n + g + h), n = 0,1,2,..., g,h>0, rj(x) = cos2x, (3.19) 
(f) n (x) = (smx) 9 (cosx) h Pji a ' (3) (cos2x), a = g-l/2, fi = h-l/2, (3.20) 

in which Pn*'®(j)) is the Jacobi polynomial (|C.8|) . It is straightforward to evaluate the 
Poisson brackets 

{H, cos 2x} = 2p sin 2x, (3.21) 
{H,{H, cos2x}} = -cos2x8H' -4(g 2 -h 2 ), H' = H + (g + hf /2, (3.22) 
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leading to the solution of the initial value problem: 



cos 



2x(t) = (cos2x(0) + 9 — ) cosfcty^] 



2H' 

-p(0) S in 2 ,(0)^£a-^. (3.23) 

Note that | cos2a;(t)| < 1 is satisfied. The corresponding quantum expressions are 

[ft,cos2x] =2(isin2xp + cos2x), (3.24) 

[ft, [ft, cos 2x] ] = cos 2x (8ft' -4) + 4 [ft, cos 2a;] + 4(a 2 - p 2 ), (3.25) 



a±{H) = 2±2v / 2ft 7 . (3.26) 
The exact operator solution reads 

e ttH cos 2x e~ ttn = (i sin 2x p + cos 2x) 



ia+(H)t _ ia-(H)t 2 _ «2 



272ft 7 2ft - 1 

+ (cos 2, (2W - 1) + a = - ^) _ (_ - _) . (3.27) 
The annihilation and creation operators are 

a' {±) /2 = a {±) 2Vm' = ± sin 2x— + cos 2x VW' + ; H . (3.28) 

dx v^ft 7 ± 1 

It is now obvious that they (a'^) are not hermitian conjugate to each other. When applied 

to the eigenvector <p n as 2£ n + (g + h) 2 = {2n + g + h) 2 , we obtain: 

a 2 -(3 2 ± 



a' { V 2 0n = -sin2x— - + (2n + g + h) cos2x0„ + ; . ■„ 

dx 2n + a + p 

4(n + a)(n + /?) 



2n + a + P 

^ a 2 - (3 2 



(3.29) 



a' (+) /2 n = sin 2x-p^ + (2n + g + h) cos 2x <p n + 



v 3 ' r 2n + a + /3 + 2 

4(n + l)(n + a + /3 + 1) 

2n + a + /3 + 2 ^ n+1 ' 



(3.30) 



The right hand sides are the results of the application. From ()3.29|) and a'( >(j) n = A(2n + 
g + h)a^(p n , the coherent state ip and ■0' are 

ij(x) = o (x) £>/2)" 7 ^TTH^ P^\cos2x) } (3.31) 

J^jJ (« + l)n(P + l)n 

^(s) = o (x) ^(A/4)" ^ 2 7 |" Pi^)(cos2x). (3.32) 
^ (a + l) n (/^ + l) n 

We are not aware if concise summation formulas exist or not. 
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3.1.3 Soliton potential, or the symmetric Rosen-Morse potential 

As is well-known — g(g + l)/cosh 2 x potential is reflectionless for integer coupling constant 
g, corresponding to the KdV soliton. It has a finite number 1 + [g\ (the greatest integer not 
equal or exceeding g) of bound states: 

7i = f (p + ig tank x)(p — igtanhx)/2, — oo < x < oo, g > 0, (3.33) 
S n = n(-n/2 + g), n = 0, 1, . . . , [g]' , rj(x) = sinhx, (3.34) 
n (x) = r n {coshx)- 9 Pi m {ismhx), (3 = -g - 1/2. (3.35) 

These eigenfunctions are real due to the parity P^\—x) = (—l) n Pn ,a \x). The Poisson 
brackets are 

{ft, sinhx} = -pcoshx, {H, {ft, sinhx}} = sinhx 2ft', H' = H-g 2 /2, (3.36) 
leading to the solution of the initial value problem: 

sin ^-2^ 



sinhx(t) = sinhx(0)cos[t v / -2^] + p(0) coshx(O) 1 V = J . (3.37) 

^-2H' Q 

It describes sinusoidal motion for bound states H' < only. But the above expression is 
valid for the unbound motion H' > 0, too. 

The corresponding quantum expressions are 

[ft, sinh x] = —i cosh x p — sinh x/2, (3.38) 
[ft, [ft, sinh x]] = - sinh x(2H' + 1/4) - [ft, sinh x] , (3.39) 
a±(H) = -l/2± V-2ft'. (3.40) 

The exact operator solution reads 



e ltn sinh x e * m = (— i cosh x p - sinh x/2) 

2 V /Z 2ft 7 

2ft' + l/4/e ia +( w )* e <«-(w)*N 
2 v / =2ft 7 i^) ~ ^H))' ( ' } 



— sinh x 

The annihilation and creation operators are 



a'W = a (±) 2 v /= 2ft 7 = =f cosh x-^- + sinh x v /= 27f. (3.42) 

ax 
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When applied to the eigenvector <f> n , we obtain as 2£ n — g = —{g — n) : 



a 



'(-), 



cosh x + (g — n) sinhx0 n = (n + (3)4> n -\, (3.43) 



dx 



, M . . dd> n , \ • i (n + l)(n + 2/3 + 1) , 

a w 0n = - coshx— h [g - n) smhx<p n = n+ i. (3.44) 

Gfe n + p + 1 

We obtain the following interesting commutation relations: 

[W,a /(±) ] = ±(v /Z Wa /(±) + a ,(±) v /Z 27f )/2, (3.45) 
[a'H a '(+)] = 2v /= 27f, a'^a'W + a'f+'a'H = 4H + 2g, (3.46) 

which look very similar to those for the 1/ sin 2 x potential ()2.65j) - (J2.b7j) . In contrast to the 
1/ sin 2 x case, the present case has only finite dimensional representation, n — 0,1, ... , [g]' , 
that is from the ground state to the highest level. There is no coherent state as the eigenvector 
of the annihilation operator 1)3.42)1 . 

3.1.4 Morse potential 

This is another well-known example of exactly solvable potential with a finite number of 
bound states 0: 

H d = (p + ifie x -ig)(p-ifie x + ig)/2, -oo < x < oo, /i,g>0, (3.47) 
£ n = n(-n/2 + g), n = 0, 1, [g]', rj(x) = e~ x , (3.48) 
(j) n {x) = e~^ eX+9X e- nx L^ 9 - 2n \2fie x ). (3.49) 

The Poisson brackets are 

{H,e~ x } =pe~ x , {H,{H,e- x }} = e~ x 2H' + /ig, H' = H-g 2 /2, (3.50) 

leading to the solution of the initial value problem: 

«-« = + $7) - [ty/=*Ri] - ,[0)^ ^g - ||. (3.51) 

It describes sinusoidal motion for bound states Ti.' Q < only. But the above expression is 
valid for the unbound motion H' > 0, too. It is easy to verify e~ x( -^ > 0. 
The corresponding quantum expressions are 

[H,e- X ] =ie- x p-e~ x /2, (3.52) 
[H, [H, e' x ] ] = -e- x (2H' + 1/4) - [H, e~*\ - fi(g + 1/2), (3.53) 



a±{H) = -1/2 ± V-2H>. (3.54) 
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The exact operator solution reads 



e itn e -x e ~itn 



(ie- x p-e~ x /2) 



i*+(K)t _ e i*-(H)t fj,(g + 1/2) 



2 V /Z W 2W + 1/4 

i , p ia + {H)t ia-(H)t. 

- ( e -(2tf + 1/4) + „( s + 1/2)) - _) . (3.55) 

The annihilation and creation operators are 

a'W = a&2V=TK = ±e~ x 4- + e^V=2W - ^£±lj . (3.56) 
When applied to the eigenvector <j> n , we obtain as 2£ ra — g 2 = — {g — n) 2 : 

a'< >4> n = -e x — + (g - n)e >„, - ~p rTT^ = ^77 n , > n-i, 3.57 

ax 2(p — n) + 1 2\g — n) + 1 

a ,(+) 0n = e + (p - n)e x <f) n - PV \ 6 n = \ \ : <W 3.58 
ax 2(g — n) — 1 2 (a — nj — 1 

3.2 'Discrete' quantum mechanical systems 

For specifying the dynamical systems belonging to the 'discrete' quantum mechanics [2113]; 
we use the name of the polynomial eigenfunctions for want of universally accepted naming. 
The factorised Hamiltonian is given by (|2.76|) . 

3.2.1 Continuous Hahn polynomial (special case) 

The factorised Hamiltonian of the continuous Hahn polynomial (special case) has a potential 
function V depending on two parameters: 

V(x) = (ai + ix)(a 2 + ix), — oo < x < oo, a 1; a 2 > 0. (3.59) 

The eigenvalues and eigenfunctions are: 

£ n = n(n + 2a 1 + 2a 2 -l)/2, n = 0,1,2,..., (3.60) 



4>o( x ) = y Ilj=i r ( a j + ix)T(aj - ix), tj(x) = x, (3.61) 
(f) n (x) = <j) (x)P n (x), P n (x) = p n (x ; a u a 2 , a u a 2 ), (3.62) 

in which p n (x ; a±, a 2 , a±, a 2 ) is a special case of the continuous Hahn polynomial (|C.11[) . This 
is a two parameter deformation of the Hermite polynomial. Thus this dynamical system is 
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a deformed oscillator. The classical solution shows this fact clearly: 

{H,x} = -^(al + x*)(al + x 2 )swhp, {H, {H, x}} = -x(2H + ( Ql + a 2 f) , (3.63) 
x(t) = x{0) cos[t v / 2Wo + (ai + a 2 ) 2 ] 

+ V(.; + «(0)»)(«l + x(0)») S inhp(0) " [ '^ M °f l (3.64) 
v V 2 ^o + (ai + a 2 y 

The corresponding quantum solution is also simple: 

[H, x] = -i(T+ - T_)/2, (3.65) 
[W, [W, x] ] = x(2H' - 1/4) + [W, x], 2H' = 2H + (tn + a 2 - 1/2) 2 , (3.66) 

e ltn x e~ ltn = [H, x] == + x — ; + v ; , 3.67 

a±(H) = 1/2 ± Vm'. (3.68) 

The annihilation and creation operators are: 

a '(±) = a^^Vm* = ±[H,x]Txa T (H) = =p(T + — T_)/2 + x(V2jj/ =f 1/2). (3.69) 

When applied to the eigenvector <fi n , we obtain as 2£ n + (a 1 + a 2 — 1/2) 2 = (n+Oi + a 2 — 1/2) 2 : 

2a' ( -V„ = - T_)0„ + 2:r(n + a, + a 2 )0„ 

= (n + ai + a 2 - l)(n + 2ai - l)(n + 2a 2 - l)0„_i, (3.70) 

2a' (+) 0« = - T_)0„ + 2x(n + ai + a 2 - 1)0„ 

(n+ l)(n + 2ai + 2a 2 - 1) 



n + a\ + a 2 
The similarity transformed operators act as 



0„+i. (3.71) 



<P (x)- l a'^<p Q (x)-P n (x) = (x(n + a, + a 2 - \ T \) =F ^(V(x)e p - V(x)*e~ p )^j P n (x).(3.72) 
The coherent state tp and ip' are 

" (2A) 

n=0 (2ax) B (2 

We do not know if these sums have concise expressions or not. 



^ (2ai)„(2a 2 )„(ai + a 2 ) n 
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3.2.2 Continuous dual Hahn polynomial 

The continuous dual Hahn polynomial has three parameters (ai, 03) and is considered as 
a two parameter deformation of the Laguerre polynomial Ln . The factorised Hamiltonian 
of the continuous dual Hahn polynomial has a potential function V: 

V \ x ) = o- /o- — r~rT~' < a; < 00, ai, a 2 , a 3 > 0. (3.75) 
2zx(2zx + 1) 

As a dynamical system this is a deformed Calogero model, or a deformed x 2 + 1/x 2 potential. 
Like the Calogero model it has a linear spectrum and the eigenfunctions are: 

S n = n/2, n= 0,1,2,..., (3.76) 



■ (x) = \ 1 J " n \ ~ 3 ~LW , V (x)=x 2 , (3.77) 



H i =1 T(aj+ix) lli=i r ( Q j - 
r(2iac) T(-2zx 

6 n (ar) = cf)o{x)P n (x 2 ) ) P n (r)) = S n (r) ; a x , a 2 , a 3 ), (3.78 



in which S n (j] ; ai, 02, 03) is the continuous dual Hahn polynomial (|C.12|) . For deriving the 
classical solution, let us note that the quantum potential ()3.75|) has acquired quantum cor- 
rections from the classical one: 



(ZIX) 



The classical motion is simple: 



nLK 2 +z 2 ) 

{H c , x 2 } = — * — ^ sinhp, (3.80) 

{H c , {H c , x 2 }} = -x 2 /4 + m\ + h H c + 62/4, h = a i > fe 2 = ^ , (3.8I) 

l<jr<3 l<i<fc<3 

x 2 (t) = (x 2 (0) - 8H 2 c0 - AhHco - b 2 ) cos[t/2] + 8ft 2 + 46^^ + 6 2 

+ -i — f- sinhp(O) sin t/2]. (3.82) 

x(0) 

The quantum version is almost the same: 

[H, x 2 } = -ix{T + - T_) - (T+ + T_)/2, (3.83) 
[ft,[ft,x 2 ]] =x 2 /4 + J R_ 1 (ft), R^(H) = -(2H 2 + (6 X - 1/2)W + 6 2 /4), (3.84) 
jtH x 2 e -m = 2i[ft,x 2 ]sin[t/2] + (x 2 + 4R_i(W)) cos[*/2] - 4i2_i(W). (3.85) 
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The annihilation and creation operators are: 

a {±) = ±[H, x 2 } + x 2 /2 + 2R- 1 (H) 

= Tix(T+ - T_) =F (T+ + T_)/2 + x 2 /2 + 2R- 1 (H). (3.86) 

When applied to the eigenvector <f> n , we obtain: 

= -n | f (n + % + a* - 1) • n _i, (3.87) 

l<j'<fc<3 

a (+) 0n = -0n+i- (3.88) 
The similarity transformed operators are: 

M^)' 1 « (±) M x ) = x2 / 2 - 47 ^ 2 - 2 ( fo i - V 2 )ft - V 2 

=F ((1/2 + zx)y (x)e p + (1/2 - ix) V(x)*e' p ) , (3.89) 

in which H = (p^x^H^x) = (V(x)e p + V(x)*e~ p - V(x) - V(x)*)/2 is the Hamiltonian 
counterpart at the polynomial level satisfying 7iP n (x 2 ) = n/2P n (x 2 ), see Appendix C. The 
coherent state is 

1>(x) = Mx) £ m 7 T s ^n(x 2 ). (3.90) 

We do not know if this sum has a concise expression or not. The commutation relations 
among TC, and a^' are more complicated than su(l, 1): 



32ft 3 + 24(6i - 1 /2)H 2 + (4(6i - 1 /2) 2 + 46 2 + l)ft + 6i6 2 - cna 2 a 3 . (3.91) 



As in the Meixner-Pollaczek case ()2.94|) . the annihilation and creation operators for the 
continuous dual Hahn polynomial factorise into the operators A and A' appearing in the 
Hamiltonian H = A^A/2: 

a H = lU, a (+) = A ] X. (3.92) 
The operator X in this case reads 

n? = i(2aj - 13 



/ FT* i(2a,--l)\ 

+ + ( x + -1) + i \. J )S_. (3.93) 

V z 8(1 + arj / 

These X and compensate the shift of the parameters (01,02,03) caused by A* and A, 

respectively. See Appendix B for more details. 
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3.2.3 Wilson polynomial 

The Wilson polynomial has four parameters (a^, a 2 , a 3 , a 4 ) and is considered as a three pa- 
rameter deformation of the Laguerre polynomial . The factorised Hamiltonian ()2.76j) of 
the Wilson polynomial has a potential function V: 

nt=i( a i + i x ) 

V(x) = 2ix ^ 2ix - ^ , <x <oo, ai, a 2 , a 3 , a 4 > 0. (3.94) 

As a dynamical system this is another deformation of the Calogero model, or a deformed 
x 2 + 1/x 2 potential. The spectrum is now quadratic in n and the eigenfunctions are: 

£ n = n ( n + Y: 4 j=1 a j -l)/2, n = 0,1,2,..., (3.95) 



M*) = \ J ~l,l\ ~ S ~L ' s L , V(x)=x 2 , (3.96) 



n* =i T ( a 3 +ix ) nl=i r ( Q j - w ) 

r(2wc) r(-2ia?) 
,{x) = 4> (x)P n (x 2 ), P n (r)) = W n (r] ; a x , a 2 , a 3 , a 4 ), (3.97) 



in which W n (j] ; O^j (22) ^3) ^4) IS the Wilson polynomial ()C.13|) . The classical motion looks 
like a cross between those of the continuous Hahn and the continuous dual Hahn potentials 
with the classical potential V c : 



{K,x 2 } = - V 3 ^ sinhp, V c (x) = n ^ x)2 (3.98) 

{H c , {H c , x 2 }} = -2x 2 (H c + ci) - R-i(H c ), (3.99) 

R^(H c ) = -2(H 2 c + c 2 H c + c 3 ), Cl = 6?/8, c 2 = 6 2 , c 3 = Ms/4, (3.100) 

-R-l(^co) \ r, ; Tl -R-l(^co) 



* 2 (*) = (* 2 (0) + ± ll 7^\ ) co 8 [V2(Wco + Cl )] 



2{H cQ + Cl )J /J 2(K + ci) 

v 'PI^ + ^» s.4*V^T7Q] , (3 , 01) 

2^(0) " ' V2(W,o + 30 



where we use the abbreviation 

& i = *V] % ' & 2 = f a ^ afc ' &3 5Z a i a *= a /- (3.102) 

l<j<4 l<j<jfc<4 l<j<fc</<4 

The quantum version has almost the same form with quantum corrections in the coeffi- 
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cients ci, C2 and C3: 

[H,x 2 ] = -ix{T + -T_) - (T+ + T_)/2, (3.103) 

[7Y,[?i,x 2 ]] = [ft, a; 2 ] + 2a; 2 (ft + c 1 ) + R^H), (3.104) 

= -2(ft 2 + c 2 ft + c 3 ), (3.105) 

Cl = &!(&!- 2)/8, c 2 = 6 2 -6 1 /2, c 3 = (6i - 2)63/4, (3.106) 



2 v / 27f 2(W + ci) 



\ -a_(W)e ia +( w >* + a+(ft)e ia -( w )* 



+ K + 2(^)J 2\/WV ' (3 ' 107) 

a± (ft) = 1/2 ± v 7 ^, 2ft' d = 2ft + 2d + 1/4. (3. 108) 

The annihilation and creation operators are: 



a'W = a^2vW = ±[ft, a: 2 ] =F :r 2 a T (ft) + 



= T ix(T + - T_) =F (T + + T_)/2 =F * 2 a T (ft) + . (3.109) 

When applied to the eigenvector n , we obtain as 2£„ + 2ci + 1/4 = (2n + 61 — l) 2 /4: 

° Vn = - (2n + 6r-2)(2n + 6l -l) ^ ' (3 - 110) 

<* ,(+ Vn = - , 9 r^TWrrrr ^- ( 3 - m ) 

(2n + 61 — l)(2n + b\) 
The similarity transformed operators act as 



4> Q {x)- 1 a! { &4> Q {x)-P n {x 2 ) (3.112) 

-R-l(^n) 
■'nil 

The coherent state ip and ■?/>' are 



= (±(£ n - S nTl )x 2 ± f R ' l(g "J =F ((1/2 + ix)V(x)e p + (1/2 - ix)V (x)* e~ p )) P n (x 2 ) . 

V W+1 — On ' 



m = W) E ( "f n + " 2 + r? t!f ' p "<" 2 >- < 3 ' 113 > 
»■(») = M*) E ( " 2 f " 2 + ? + ° 4 '( 2 '° p "(" 2 )- < 3 ' 114 > 



n=0 



n! ni< i<fe <4( a i + a k)n 
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4 Summary and Comments 

Unified theory of annihilation-creation operators a& is developed for various exactly solv- 
able quantum mechanical systems possessing the 'sinusoidal coordinate'. It applies to most 
of the degree one solvable quantum mechanical systems as well as the solvable 'discrete' 
quantum mechanical systems, which are also shape-invariant j2]. The eigenfunctions of the 
latter are described by the Askey-scheme of hypergeometric orthogonal polynomials j3j. The 
method provides an independent algebraic solution of these quantum systems. The energy 
spectrum is obtained a la Heisenberg and Pauli from the Heisenberg operator solution for 
the 'sinusoidal coordinate' i], e rje~ and the entire eigenfunctions are explicitly obtained 
as {(a( + )) n 0o}, n — 0, 1, . . ., in which 0o is determined by aS~'(j)Q = 0. Various examples 
are worked out in section |2] and El It also applies to theories with a finite number of bound 
states. It should be stressed that these annihilation-creation operators are natural ones 
containing the differential (difference) operators, in contradistinction to those annihilation- 
creation operators introduced in the algebraic theory of coherent states jH]- By a similarity 
transformation in terms of the ground state wavefunction <f) , the Heisenberg operator so- 
lution gives the structure relation for the corresponding orthogonal polynomials j3] and the 
annihilation-creation operators provide their shift down-up operators. Another characteristic 
feature is the uniqueness. Except for the overall factor, which is intrinsically undetermined, 
the action a^(j) n is completely determined by the Hamiltonian of the system. This means 
that the relative weights of the terms in a^<p n are governed by the energy spectrum. We 
have shown in some detail that this type of algebraic exact solvability is valid at both clas- 
sical and quantum levels. This is in good contrast with shape-invariance, which is a strictly 
quantum notion. The necessary and sufficient condition for the existence of the 'sinusoidal 
coordinate' is worked out for the ordinary quantum mechanical systems in Appendix A. It 
is a good challenge to derive a corresponding result for the 'discrete' quantum mechanical 
systems. 

Generalisation of the present formalism to multi-particle systems is highly desirable. 
Simplest multi-particle systems possessing the 'sinusoidal coordinate' and the corresponding 
Heisenberg operator solution is the Calogero systems based on any root system ^T] . In fact 
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a more general Hamiltonian 

n n d 

n = -J2(p" + ^) + nx), V(x) = -2V(x), (4.1) 

3=1 j=l 3 

n 

[H,[H, V ]]=A( v -n), v = J2 x l ( 4 - 2 ) 

3=1 

of harmonic oscillators modified by a generic homogeneous degree -2 potential has the same 
property. The corresponding eigenfunctions are the Laguerre polynomials again [To"! ITT] . 
As is well known the annihilation-creation operators of the harmonic oscillator have a quite 
wide applicability in many branches of physics. We wonder if the newly found annihilation- 
creation operators for the other solvable quantum mechanical systems might find an equally 
wide range of applications. 
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Appendix A: Determination of the potentials having the 
'sinusoidal coordinate' 

We have seen that the existence of the 'sinusoidal coordinate' or the exact Heisenberg oper- 
ator solution (|2.14|) leads to the unified definition of the annihilation-creation operators. All 
the examples discussed in the text share the common property of 'shape-invariance', thanks 
to which the corresponding quantum systems are exactly solvable. Here in Appendix A 
we analyse, within the context of ordinary quantum mechanics, the necessary and sufficient 
condition for the existence of the 'sinusoidal coordinate' and show that such systems consti- 
tute a sub-group of known 'shape- invariant' quantum mechanics. For the 'discrete' quantum 
mechanical systems, writing down corresponding conditions is easy. It would be a good 
challenge to provide a complete list of 'discrete' quantum mechanical systems admitting the 
'sinusoidal coordinate' or the exact Heisenberg operator solution. 

For a given pair (j](x),7i) of a coordinate function rj(x) and a Hamiltonian 7i to satisfy 
the exact Heisenberg operator solution (|2.14|) is equivalent to the condition (|2.12|) that the 
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multiple commutators of H with r\ form a closed algebra at level two 



{&dH) 2 ri(x) = [n,[n,v(x)}} = vRo{n) + [n,v(x)}R 1 (n) + R- 1 (n). (a.i) 

Here the coefficients Ro(7i), R\{7i) and R-\{7i) are polynomials in the Hamiltonian 7i. 
It should be stressed that this condition is purely algebraic and the knowledge that the 
eigenfunctions have the general structure (|2.3|) is irrelevant. The latter ()2.3j) is a consequence 
of the condition (|A.1J) . For the ordinary quantum mechanical system with potential V(x) 

the commutator between TC and rj reads 

^^-V'^-Iv", (A.3) 
(adW) a ^ = ^ + ^ + i»T + ^ / , (A.4) 

in which primes denote differentiation with respect to x. From (|A.4J) we see that the l.h.s. of 
flA.I|) contains the derivative operator (the momentum operator) at most quadratic degree. 
So must be the r.h.s. since the momentum operator can come in as a part of 7i (as p 2 /2) or 
as [H, rj\, see flA.3|) . Then we can parametrise 

R (H)=r^H + r^\ R 1 (H) = r u R^{U) = r^H + r% (A.5) 

in which rj are all constants. Then the coefficients of the operators ^ and the function 
part of (|A.1|) give the conditions: 

^ — r ^ (A 7) 

doc doc 

1 A dr L dV_ _ 1 A , f (i) , (iK y , (o) (o) , A , 

The first condition (jA.6|) simply means that 77 is either a trigonometric or a hyperbolic 
function of x which gives an exponential function or a quadratic and linear polynomial in x 
in the degenerate limits. By comparing (jA.6|) and ()A.7|) we obtain 



rj 1} = 2r x . (A.9) 
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Then (jA.8|) reduces to 

drjdV d 2 r] / 1 \ {0) (0) 
which integrates easily when multiplied by drj/dx: 

Here c is the constant of integration. Thus we have determined the possible form of the 
potential V in terms of the 'sinusoidal coordinate' 7](x) and its first derivative drj/dx, with five 
parameters r%, r^, r^j, and c in (jA.lOj) and two more possible constants of integration 
of (lA~6j) : 

1 (o) 

VV) = 7^^ + r% + c)-%. (A.ll) 



The actual number of essentially free parameters is much less, since the origin of the quadratic 
potential, or the location of the singularity, etc, could be freely adjusted by introducing new 
variable x new = ax + (3 (H new = H/a 2 ). The condition that the Hamiltonian must be 
bounded from below imposes some constraints on the parameters. The overall additive 
constant is fixed uniquely when the ground state energy is required to be vanishing £ = 0. 

It is rather straightforward to determine all the potentials possessing the 'sinusoidal 
coordinate' and thus algebraically exactly solvable. They all belong to the known group of 
shape-invariant potentials. Except for the trivial case V = 0, we have 

1. Rational case, r\ = 0. The generic solution of (jA.6|) is 



rj(x) = — -r^lx 2 + C\X + c 2 , (A. 12) 

with ci and c 2 being the constant of integration. Two special cases are of interest: 
rj(x) = x gives the harmonic oscillator and rj(x) = x 2 leads to the x 2 + 1/x 2 potential 
discussed in section 13.1.11 



2. Trigonometric case, r\ > 0. The generic solution of (jA.6|) is 

r (i) _ 
7](x) = — — — + Ci cos ^fr{x + c 2 sin y/r^x, (A. 13) 

Z7"l 

with C\ and c 2 being real constants of integration due to the reality (hermiticity) of 
r). By rescaling and shift of the coordinate x, it reduces to the Poschl- Teller potential 
discussed in section 13.1.21 The 1 / sin 2 x potential in section 12.1.11 and the symmetric 
top are obtained as degenerate cases. 
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3. Hyperbolic and exponential cases, r% < 0. The generic solution of (jA.6|) is 

r (i) 

r](x) = h ci cosh \J —r\ x + C2 sinh a/ - r i x > (A. 14) 

in which the constants of integration c\ and C2 could be vanishing or equal c\ = ±C2. 
The generic case leads to the hyperbolic Poschl- Teller potential and the degenerate 
cases contain the soliton potential in section 13.1.31 and hyperbolic symmetric tops 
and the Morse potential in section I3.1.4[ etc. We could not discuss all due to space 
limitation. 

In all these examples, the prepotential W has also a simple expression in terms of 77 and 
drj/dx: 

ax 

Here the prepotential W is related to the ground state wave function <p and thus to the 
potential V as 

-K(§) 2 ^ 

and it plays an important role in supersymmetric (shape-invariant) quantum mechanics 

muni 

It should be stressed that not all shape-invariant and exactly solvable potentials admit 
the 'sinusoidal coordinate'. Such examples are the Kepler problems in rational, spherical 
and hyperbolic coordinates and the Rosen-Morse potential, respectively: 



"M = ±(-U4q^ + i), (A.16) 

2 V x ar g z 1 



V(x) = \ (-2n cot x + rt^- + ^-9 2 ), (A. 17) 
2 V sin x g z / 

= V-2/icothx + 9 -^fi- + ^+g 2 ), (A. 18) 
2 V sinn x g / 

7(ar) = i^tanhs - + ^ + /)■ (A. 19) 
2 V cosh x g z / 

Their wavefunctions do not have the general structure (|2.3jl . either. 
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Appendix B: Interpretation in terms of Shape Invari- 
ance 



As shown in section 121 the annihilation-creation operators are completely determined once 
the closed relationship (j2.12j) among rj, [H,rj\ and [H, [H,rj\] is obtained. Although it plays 
no active role in the determination of the annihilation-creation operators, shape invariance 
is the common property underlying all these exactly solvable Hamiltonians discussed in this 
paper. Therefore it is interesting as well as illuminating to understand the mechanism of the 
annihilation-creation operators within the framework of shape invariance. For this purpose 
we concentrate on the annihilation-creation operators of the Meixner-Pollaczek polynomials 
(J2.94|) and of the continuous dual Hahn polynomials (|3.92|) . which factorise into the operators 
A and A^ constituting the shape invariant Hamiltonian Ti. = A' A/2. Another motivation 
of this Appendix is to provide a bridge between the physics of 'discrete' quantum mechanics 
[2] and the analysis of Askey-scheme of hypergeometric polynomials |BJ . The latter focuses 
on the polynomial part of the eigenfunctions, whose orthogonal measure is provided by the 
ground state wavefunction (|2.4|) . 

Let us start with recapitulating the rudimentary facts of the shape-invariant 'discrete' 
quantum mechanics as developed in |2J. Knowledgeable readers may jump to the main 
results flR2^ - (lR29|) . but some intermediate results (li3~7l) and flB7T3l) flB7T7|) would also be 
interesting in connection with the 'sinusoidal coordinate' rj(x). 

A shape invariant quantum mechanical system consists of a series of isospectral Hamil- 
tonians {H(X)} parametrised by (a set of ) parameters A = (Ai, A2, • • • ): 

H{\) = ^l(A)U(A)/2, 0„(x;A) = o (x;A)P„,(r / (a;);A), E n {\), etc. 

Shape invariance is tersely expressed as 

^(A)^(A)t = .4(A + <5)U(A + 8) + 2£ X (A), (B.l) 

where 8 is a shift of the parameter. 1 The operator .4(A) maps the eigenvectors of TC(X) to 

those of 7^(A + 8) and the other operator *A(A)^ acts in the opposite direction. In each case 

1 In the case of the Askey- Wilson polynomials, this is modified to A(X)A(\y — q 2S A(\* q s )^A(\* q s ) + 
2£i(A), where 5' is a constant and A * q s = (\iq Sl , \2q S2 , • ■ •)■ 
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studied in this paper, the parameter A and the shift 8 are 



Meixner-Pollaczek : 


A = 


a. 


8 


= 1/2, 




(B.2) 


continuous Hahn : 


A = 


Oi- 


a-i, 


, s = 


;i/2,i/2), 


(B.3) 


continuous dual Hahn : 


A = 


(01, 


a 2 . 


a 3 ), S 


= (1/2,1/2,1/2), 


(B.4) 


Wilson : 


A = 


(a l5 


a>2, 


a%, a^j, 


S = (1/2, 1/2, 1/2, 1/2), 


(B.5) 


Askey- Wilson : 


A = 


(ai, 


02, 


as, 04), 


6 = (1/2, 1/2, 1/2, 1/2), 


5' = -1/2. (B.6) 



The ground state 4>o(x) and the orthogonal polynomial P n (r)(x)) are given in (|2.83|) - (|2.84j) . 
$M-$M, dHZD-(HIHD, and CTH-dZM). These ground states satisfy 2 



0o(z- i/2; A + 5) = 




X)(p(x-i/2)(f)o(x;\), 


) is given by 






Meixner-Pollaczek 


■ <p(x) 


= 1, 


continuous Hahn 


■ f( x ) 


= 1, 


continuous dual Hahn 


■ f( x ) 


= 2x, 


Wilson 


■ f( x ) 


= 2x, 


Askey- Wilson 


■ f( x ) 


= —2 sin x = i(z — z~ x 



(B.7) 

(B.8) 
(B.9) 
(B.10) 

(B.H) 
(B.12) 



Let us consider S±(X), T±(A), .4(A) given in ()2.77|) - ()2.79|) . By using the property ()B.7|) . 
we have 

M x ; A + S)- 1 S ± (X) o (x ; A) = ^(x)" 1 e ±p ' 2 , (B.13) 

0o(x;A)-^ ± (A)t o (x; X + S) = l^'^l^] , (B-14) 

[V(x ; A)* e p < 2 ip{x) . 

(In the case of the Askey- Wilson polynomials, the following replacement is needed: A + d =>- 
A * q 5 , e ±p/2 q ±D/2 , V(x ; A) V(z ; A). Hereafter we will omit similar remarks.) From 



2 For the Askey- Wilson polynomials, this relation reads 4>q{x — ij/2;X * q s ) = \/V{z ; A) ip(x 
ij/2)(j) (x ; A), where 7 = logg. 
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this, we obtain 

F(X) d = o (x ;X + d)~ 1 A(X) o (a;;A) =^(x)- 1 (e p / 2 -e" p / 2 ), (B.15) 
B(X) = <j) (x ; X)- 1 A(\) j o (a;;A + (5) = -i(V(x ; X) e p/2 - V(x ; \)* e- p/2 )<p(x) , (B.16) 

T±(X) = </>o ; A) -1 7±(A) o (x;A) = r[ X; ^f (B.17) 

Therefore the similarity transformed Hamiltonian is 

W(A) = o (x ; A)- 1 W(A) o (x ; A) = B(\)F(X)/2 

= (f + (A) + f_(A) - V(x ; A) - V(x ; A)*)/2 . (B.18) 

which acts on P n {j]{x) ; A) as H(X)P n (r)(x) ; A) = £ n (A)P n (77(a;) ; A). 

The forward shift operator P(A) and backward shift operator -B(A) act on P n (?7 ; A) as 

P(A)P n (r? ; A) = f n (X)P n ^( V ; A + 8), (B.19) 

B(X)P n ( V ;X + S) = b n (X)P n+1 ( V ; A), (B.20) 

where / n (A) and b n (X) are constants satisfying the relation / n (A)fe n _i(A)/2 = £ n (X): 

Meixner-Pollaczek : / n (A) = 2, 6 n (A) =n + l, (B.21) 

continuous Hahn : / n (A) = n + 2a\ + 2a 2 — 1, b n (X) — n + 1, (B.22) 

continuous dual Hahn : /„(A) = — n, 6„(A) = —1, (B.23) 



Wilson : /„(A) = -n(n + E •=!%-!), 6 n (A) = -1, (B.24) 

Askey-Wilson : / n (A) = -<f /2 (<T" - 1)(1 - a^a^- 1 ), b n (X) = -q-^/ 2 . 

(B.25) 

For the Meixner-Pollaczek and the continuous dual Hahn cases, we have seen that 
the annihilation-creation operators are factorised oc ATU and oW oc „4 f X, (l2~Mll . 
(f3~32l . By using (pl3|) - (lB~T4j) and (pTffl . O (^ ; A + S)~ 1 X(X)(f)o(x ; A) and o (x ; A)" 1 
A" (A) Vo(^ ; A + 6) can be written down explicitly They act on P n {i] ; A) as for the Meixner- 
Pollaczek case: 

o (x ; A + 8)- 1 X(X) M% ! A) • P n (v \ A) = 2P n ( V ; A + <5), (B.26) 
0o(x ; A)" 1 X(A) f o (x ;X + d)-P n (j 1 ;X + d) = (n + 2a)P n ( V ; A), (B.27) 
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and for the continuous dual Hahn case: 



; A + S)- 1 X(X) (f> {x ; A) • P n (n ; A) = P n (r) ;X + S), (B.28) 

o (x;A)- 1 X(A) t <f> {x;X+d)-P n {r];X + 5) = J] (n + aj + a k ) ■ P n (rj ; A). (B.29) 

i<i<fc<3 

Therefore (X) compensates the parameter shift caused by A (A* 1 ), so that the effect of 
a 1 -™-* (a^) is to give the polynomial with the same parameter A of degree one lower (higher). 
This result is new. 

Let us close this Appendix with a remark on the formal definition of the annihilation- 
creation operators used within the framework of shape-invariant quantum mechanics 
E]. A unitary operator U (W) is defined as a map between two orthonormal bases with 
neighbouring parameters, {4> n (x; A)} and A + <5)}: 

UMx;\) = Mx;X + 6), tf(j> n (x ;X + S) — 4> n (x ; A). (B.30) 

This allows to introduce new annihilation-creation operators in a factorised form 

a = Z/U(A), a) = A(XyU, (B.31) 

which satisfy TC = = A{X)^A{X)/2. The operator U is rather formal and it cannot 

be expressed as a differential or a difference operator. This operator U can be considered as 
unitarisation of the natural factorisation operator X discussed above. 

Appendix C: Some definitions related to the hypergeo- 
metric and g-hypergeometric functions 

For reader's convenience we collect several definitions related to the (g-)hypergeometric 
functions [6J. 

o Pochhammer symbol (a) n : 

n 

(a) n d = JJ(a + k - 1) = a(a + 1) • • • (a + n - 1) = T(a + n)/T(a). (C.l) 

k=l 

o g-Pochhammer symbol (a ; q) n : 

n 

(a ; q) n = " a ^ = (1 - a)(l - ag) ■ ■ ■ (1 - aq^ 1 ). (C.2) 

k=l 
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o hypergeometric series r F s : 



h, ■■■ ,b s 



dcf 



E 

n=0 



(til, , <^r)n -2- 

(61, • • • ,6 s ) n ra! ' 



where (<n, • • • , a r ) n = f rp=i( a j)™ = ( a 0n • • • K)n- 
o g-hypergeometric series (the basic hypergeometric series) r (fi s : 



a f ai > 
\bi, 



,b s 



q;z 



(Ql; " " " ; a r j ?)n 

^ (61, ••• ,6 a ;g)n 



s— r)n (1+s— r)n(n— 1)/2 



where (a x , • • • , a r ; q) n = rij=i( a i 5 <?)« = («i 5 ?)n • • • (<ir 5 <?)n- 
o Bessel function J a (;z) : 



,dcf (^/2) a 



T(a+l) u "Va + l 

o Hermite polynomial H n (x) : 

H n (x) d ^ (2x)" 2 F ( 
o Laguerre polynomial L^\x) : 



Va + 1 4/ 



'-n/2, -(n-l)/2 



x 2 ) 



L( ) def (a + 1) 



~~ i i 

n! Va + 1 



o Jacobi polynomial P, 



(«./*), 



p(^) (;r) 

which satisfies P^' a \x) = (-l) n P^' p ' \-x). 
o Gegenbauer polynomial C$?\x) : 



dcf (a + l) n f—n, n + a + (3 + I 



\n p( Q ./3)/ 



/— n, n + a + 



1 — x 



)■ 



.7: 



def 



(2A) n 



(A + 1/2) B n 
o Meixner-Pollaczek polynomial Pif 1 (x',(f>) '■ 



p(A-l/2,A-l/2), x 



a + 



1-e 



-2i<^ 



o continuous Hahn polynomial p n (x ; ai, 02, a^, a 2 ) : 

/ / / x def .„ (ai + a' 1 ) n (ai + a' 2 ) n 
p n (x;ai,a 2 ,a 1 ,a 2 ) = % : 



x 3P2I 



-n, n + ai + a 2 + a' x + a' 2 — 1, ai + ix 
a\ + a' l5 ai + a' 2 
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which is symmetric under a 4 <-> a 2 and a[ <-> a' 2 separately, 
o continuous dual Hahn polynomial S^x 2 ; ai, a 2 , a 3 ) : 

-n, ai + ix, ai — ix 
ai + a 2 , ai + a 3 

which is symmetric under the permutations of (ai,a 2 ,a 3 ). 

o Wilson polynomial W n (x 2 ; a l5 a 2 , a 3 , a 4 ) : 



S , n (x 2 ;ai,a 2 ,a 3 ) = f (a 4 + a 2 ) n (a 1 + a 3 ) n 3 F 2 f n , a i + lx i a i %x j\ (C.12) 

V ai + a9. ai + as / 



W^a; 2 ; ai, a 2 , a 3 , a 4 ) = (ai + a 2 ) n (a 4 + a 3 ) n (ai + a 4 ) 



x4F3 /-n,n + E^-l,ai+^,ai-- ^ (ai3) 
V ai + ao, ai + as, ai + a 4 / 



ai + a 2 , ai + a 3 , a\ + a 4 
which is symmetric under the permutations of (a±, a 2 , a 3 , a 4 ). 
o Askey- Wilson Hahn polynomial p n (cos:r ; ai, a 2 , a 3 , a 4 |g) : 

p n (cosx ; ai, a 2 , a 3 , a 4 |g) = f a^ n (aia 2 , aia 3 , aia 4 ; q) n 



• " • '(■><l:>,<l.:<l" '. 

a 4 a 2 , a^a 3 , aia 4 
which is symmetric under the permutations of (a 4 , a 2 , a 3 , a 4 ). 



/g n , aia 2 a 3 a 4 g n x , aie", a x e ix \ 
x 403 (C.14) 
V aiao, aias, aia 4 / 
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